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SIMULATION OF FOREST INSECT PEST. MANAGEMENT 


D. Gordon Mott 


INTRODUCTION 


Interaction with a computer in a time-sharing environment permits 
the user to perform calculations in a variety of unique ways. Among 
these are simulations of the dynamic behavior of systems of concern 
to the user, together with the opportunity to perform control or 
regulatory actions during the Sauna tian | In jargon, the user is 
permitted to play "what if?" This ene permits -experimentation 
with a numerical representation of a real system rather than with the 
system itself. As a result, experiments can be conducted on the computer 
which would be impossibly expensive with the real system, or take 


generations to complete (as in the case of most natural resource systems), 


or which would be extraordinarily hazardous. 
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There is another valuable aspect to simulation programs -- they 
permit the programmer to sp iy view of the real system and 
transmit it to another person in a form which is far more effective 
than either the written or spoken word -- the material which is 
communicated can be used and manipulated by the second person, and 
thus absorbed through direct experience. This is the principal intent 
of the present program -- to communicate to those concerned with 
reguiating natural pest populations, some of the structure which mst 
be provided in making management decisions. A knowledge of the dynamics 
of the managed pest population, the field effectiveness of the control 
material, the effect of the pest on the production of host values, and 
the economic framework within which decisions must be made, are some 
of the aspects of decision making that are to be considered. It is 
generally true that for a variety of reasons, forest pest management 
is not based upon afi adequate knowledge of any of these, and forest 
pest research is not adequately funded to produce the required knowledge. 
In a very real sense, the author is playing "what if?" -- he is trying 
to communicate to the present reader what could be done if there was 
sufficient knowledge available. In order to do so he has brought 
together a number of elementary considerations with the intention of 


illustrating the existence of an opportunity to elevate the standard 


of practice in the field. 
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This interactive computer program is intended to be a 
communications device. What is tol be communicated will become 
clear below. At this point it is very important that the 
reader/user realize what is not to be communicated. The program 
does not produce simulations for any actual events. No result is 
represented as being realistic. There are several places in the 
text where a relationship is suggested .between the results or 
structure of the program and certain real situations. In no such 
case is it being suggested that there is detailed structural agreement 
with the real case, nor that the program as a whole 3 are any 
real case. In other instances, ect are made which reveal the 
nature of the discrepancy between this simulated program and certain 
known features of 2 real situation. In brief, whatever the merits 


of this program, it is not designed to obtain solutions to real 


problems. A 
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Description of Program 


General 
The simulation contains four major elements: 
(1) a simulator for the dynamics of the pest population. 
(2) a relationship between pest population density and 
intensity of damage. 
(3) a control "device". 


(4) a cost-benefit analysis. 


There are three population dynamics models, four damage functions, 
and two control procedures. Each is discussed below in detail. It 
is poSsible to choose 2h different combinations of population dynamics 


models, damage functions, and control procedures. 


The parameters of each population dynamics model are fixed. 
However, some of those of the damage functions and the control devices 
are variable, so that the user can actually explore a very much larger 


number of situations than the 2h mentioned above. 


The program is written in 360-Fortran IV and uses a plotting 
subroutine developed at the Yale Computer Center for typewriter 


terminal plotting. A 132 character terminal is necessary. 
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Population dynamics simulation 


There is a choice of three different population dynamics 
simulators -- one which produces an oscillating density over time, 


and two which produce stable densities. 


The oscillating model consists of the difference equation: 


ig 


log Ne ay = Co + Cy log Ny + Cy log Nioy + C3 log Neo 
where: Nj = Pest density in generation j. 
C; = Constant 


This model produces an oscillation in density from about 1 x 10° 
to about 6 x 10°, If initial conditions outside the normal population 
cycle are Beyadaare transient conditions may exceed this range. The 
model was derived from an analysis of actual pest density data in which 
density in previous years affects the rate of change in density in the 
current generation, because there are lag effects in the system. One 
of these was speculated to consist of deterioration in the nutritional 
quality of current foliage when defoliation has occurred in previous 
years. Others probably consist of the normal lag in parasite, predator, 


and disease effects in such systems. 
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The model oscillates with a period of about eight years. A 


typical pattern is depicted in Fig. 1.1/ 





Af Hint: The system posseses an equilibrium point at 


log Ne = aut where 
Pag =3 
Cy = +2.615 
Cy = +1.454 
C, = -0.978 
fee = -0.0173 


-Figure 1.--Logarithm of pest density over time as generated by the 


oscillating population dynamics model. 

There are two different stable models, one based upon an 
exponential density-dependent rate of change, and one based upon 
a sigmoid (logistic)model for the same process. In both cases it 
is assumed that: 


N R..N 


ew We eet 


where R represents the net rate of population change in generation t. 
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Figure 1.--Logarithm of pest density over time as generated 


population dynamics model. 
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In the first case (see Fig. 2): 


R= Maa Leen.) 
where: C= Maximum R = 2.5 
a= .9163 x 107-5 


**designates exponentiation 


a a a 


These values of C and a have been chosen in order to give a stable 


density at 100,000. The equilibrium point in this system (where R = 1) 


is at N= 1nC/a. The time trajectory through which density moves — 


to the equilibrium value, given an initial condition elsewhere, 
depends upon other properties of the relationship. However, this 
function has been chosen to represent a case in which R, the 
capacity of the pest to increase, increases (to an upper bound) 

as pest density is reduced. That is, a higher percentage control 
will be needed to maintain density at a value lower than 100,000, 
as density is reduced. ! 

Figure 2.--Density dependent population change function used in the 


exponential population dyanmics model 


R = Ce~ANt 
Gemn7s ’ 
a me oles x 107° 
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the exponential population dynamics model 
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8. 
The second stable relationship (Fig. 3) is based upon the logistic 


function in which: 





C 
R= ———_. 
1+e**(atb/N.) 
where: CaS. 0 
a= 2.197 
b = -81093.0 





Figure 3.--Density dependent population change function used in the 


logistic population dynamics model. 
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1+e** (at+b/Nz) 
Co2,5.0 
aoe 22197 
b = -81093.0 





This system again has a stable equilibrium point at N = 100,000. 
It differs from the previous system in that R behaves differently as 
density changes, and thus the same control method will be found to be 


optimal at a different density than in the previous system. 


3 


Lliw bortiom Loxjnoo ompe oft aunt bas ,zogtad> %: 


weteye euoiven¢ oft ni net? \okenob tnate2tib s ta i 












> ee as 


7 
a4 


i, 
: Ae ery 
boend 2f (& .gf8) qifienotsater old ie diacetnie 
= a i 
oe 
pioidw ab ae 
<x 


2 _ 
( Viides)**eef . 
a 


_ 


0.2*2 ” 


TOl.S =e 


0.20088- = d 
. 


oe — 















aoe! 
oq tnobnegsh ytiened--.6 


om eoimanryb aottaluqog: a 


2 _ 
——— —— erent 


(Video) **ere 


: 
3 


> = 
0.2 a a 
f01.5 =4 


O.c@018- = d 





= > “oa 

b et = rT yf ° i ? ¢ ° bi y 
Iiiog MTILL supe o1GeIe & end ciggs woteye ¢ 
: = no 

sit tag nk toseye ewolvetq of? mot? 270%; 








8a. 


, eageess |: eee OF ~ 200,000 
Pest Density 
Figure 3.-~Density dependent population change function used in 


the logistic population dynamics model. 
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The time trajectories through which population density moves 
towards the equilibrium value for each of the two stable models 


are shown in Figure 4. 


Figure 4.--Time trajectories of pest density (initial number = 1000) 
according to the exponential and logistic, stable population dynamics 


models. 


Both of these systems approximate some general population change 
-characteristics of real systems in that they exhibit a maximum 

rate of change, and an asymptotic approach to a minimal rate. 

The assumption that the R function depends upon density is an 
elementary realistic one in many cases, although somewhat controversial 
as a general case. In fact, the assumptions throughout this program 
that the relationships are deterministic, are somewhat naive. and 


elementary. 
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Figure 4.--Time trajectories of pest density (initial number = 1000) 


according to the exponential and logistic, stable population dynamics 


models. 
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10. 
Pest Density - Damage Functions 
One of four different relationships between percentage host 


damage and pest density can be chosen: 


(1) D = 100 (1 - 10**-aN) 

(2) D = 100 (1 - 10**-a log N) 
(3) D = 100/(1 + 10**a + bN) 

(4) D = 100/(1 + 10**a + b log N) 


Where D = percentage of complete damage, and N = pest density. 
‘(Functions 2 and 4 can be stated in simpler terms, they are presented 
this way to show their relationship to functions 1 and 3.) The forms 


of these functions are shown in Figs. 5 and 6. 


Figure 5.--Damage 

Damage functions 1 and 2. In each case the density at which 
90 percent damage occurs is specified by the user, and the value 
of a calculated in the program. | 

Function 1: % Damage = 100. [1 - 10**(-a x density)] 


° 


Function 2: % Damage = 100. [1 - 10**(-a x log density) ] 
The pest density scale is the same as that in Figure 6 for 


comparative purposes. 





Figure 6.--Damage functions 3 and 4. In each case the densities at 
Syhich 10 percent and 90 percent damage occurs are specified by the user. 
Function 3: % Damage = 100/[1 + 10**(a + b x density) ] 

Fimetion 42° Damage = 100/[} + 10**(a + b x log density) ] 


The pest density scale is the same as that in Figure 5 for 
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comparative purposes. ; eee eee 
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Figure 5.--Damage functions 1 and 2. 


In each case the 
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density at which 


90 percent damage occurs is specified by the user, and the value of 


a is calculated in the program. 


Function 1: % Damage 


Function 2: 4% Damage. 


i} 


-100.[1 - 10**(-ax density) ] 


100.[1 - 10**(-ax log density) ] 


The pest density scale is the same as that in Figure 6 for comparative 


purposes. 
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Figure 6.--Damage functions 3 and 4. In each case the densities at 


which 10 percent and 90 percent damage occurs are specified by the 


user. 


Function 3: % Damage = 100/[1 + 10**(a + b x density) ] 


i} 


Function 4: % Damage 100/[1 + 10**(a + b x log density) ] 


The pest density scale is the same as that in Figure 5 for comparative 


purposes. 
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In each case, the parameters a and b are selected by the user 


(see below: "Using the program'"'). 


Control Procedures 

Two different kinds of control are provided: an insecticide, 
and sterile male release. 

The insecticide is assumed to possess a typical sigmoid relationship 
between the percentage mortality procuded, and the logarithm of 
dosage applied. A logistic function is used to describe the 
relationship (Fig. 7). Its parameters are derived from specification 
of the LD50 and LD95 of the material, in response to a query from 


the program. 


Figure 7.--Dosage-response function for effect of insecticide. The 


curves exhibit the same LD50 of 2, and different LD95's of 10 and 20. 
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Figure 7.--Dosage-response function for effect of insecticide. 
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The sterile-male control procedure is assumed to have an effect 

directly proportional to the proportion of sterile males in the male 
population. The male population is assumed to be one-half of the 


total population. A "dose-response" relationship for this approach 


is shown in Fig. 8. 





Figure “gana ee na relationship for sterile male release. 
Percentage "mortality" obtained from given values of sterile male/ 
population ratios (S/D), assuming fully competitive sterile males, 
“and a .50 sex ratio. 


Percentage mortality = 100 S/[{(D/2),+ S] 


where, 
S = number of sterile males 
D = pest density 


In each case, the effects of the control procedure are assumed 

to be as follows: 
Negy = (1 - PM/100) Ni, 
where N,,, = actual number in next generation. 


Nty] = predicted number in next generation, 


PM = percentage mortality from control. 
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Figure 8.--"Dosage-response" relationship for sterile male release. 
Percentage "mortality" obtained from given values of sterile 
male/population ratios (S/D), assuming fully competitive sterile 
males, and a .50 sex ratio. 

_ Percentage mortality = 100 S/[(D/2) + S] 


where, 
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13% 
Thus, no provision is made for such known effects as increased 
natural survival in the remainder of a generation to compensate 
for control-induced mortality within the generation, as a result 
of density-dependent survival phenomena in the system. Instead, 
this assumption yields the result that is usually naively expected 
in control operations -- a degree of reduction in pest density equal 
to that obtained by control. 
The user is permitted to choose whether he will always apply 
the same level of control, or will select a unique level in each 
Pre nes ation’ In the case of sterile male control, he is permitted 
to specify a percentage mortality (the program then calculates 
dosage as the number of etetite males released), or to specify 
the number of sterile males to be released (the program calculates 


percentage mortality). 
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14, 
Cost-Benefit Analysis 

The program calculates the benefits derived from control 
activities on two different bases, termed short term and long 
term. "Short term'' gains and benefits (gain minus cost) are 
calculated from current control costs, and the difference between 
the value of actual damage, and that to be expected from the current 
uncontrolled predicted pest population. "Long term" gains and 
benefits are calculated from the current control costs (as above) 
‘and the difference between the value of actual damage, and that 
to be expected if no control activities were ever carried out. 
These two a cibtiéns are needed, because, for example, certain 
control activities will result in higher average damage over the 
long term, even though they result in positive gains af the year 
in which they are conducted. This phenomenon can be seen most 
easily in the case of the oscillating model. On the other hand, 
certain control actions will have the opposite result -- over the 
long term, average damage is reduced even where there is a net 
loss in the year of the control activity. 

In detail, cost, gains, and benefits are calculated as follows. 
prenever a control action is taken, it is assumed that a fixed 
cost (FIX); and a variable cost, the cost per unit of the control 
‘material (VAR) multiplied by the number of units used (DOSE), are 
incurred. Thus, 

COST = FIX + VAR.DOSE 
When no control action is taken (percentage mortality = 0) no cost 
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is incurred. : 
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The gain in value as a result of control action (VGAIN), is 
calculated by subtzacting the value. of the percentage damage actually 
obtained from that expected if no control action were taken. Short 
term and long term gains are each calculated as described above, using 
the two different pest densities for calculations of expected damage. 

The concept of the value of damage is purposely left very 
general. The user determines what value he will place upon 100 
percent damage -- it can be the value of a year's growth reduction, 
the value of an aesthetic impact, the value of a crop destroyed, 
“or any combination of these. Within the concept of this program, 
this does not matter. Something that does matter, however, is 
that no provision is made for the effects of damage to compound. 
That is, in real forest systems, one expects that in some cases 
severe defoliation, say, in two or three consecutive years in a 
decade will have a more severe impact than the same severe damage 


separated by years of light damage in which growth recovery takes 


place. 
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Similar remarks might be made concerning the costs of a unit 


of control material. It is not intended that this should be simply 


| 


the cost of treating with say, one pound of material per acre. The 


user can decide whether more general environmental costs are also to 


‘ 2 : : : oc 
be included. 2/ However, in this case again, no provision has been 


A A 


2/ Consideration was given to including constraints on the amount 


of any toxic material used. However, while this sort of restriction 


must be involved in any real case, the purposes of this program do not 


require it. 


made for say, compound effects from the accumulation of toxic material 


in the environment. 

Finally, no provision has been made for compounding interest on 
either the costs of control or the value of gains. This would also 
be required in any real analysis, as would provision for rising trends 
in costs and values in an economic projection. 


Using the Program 


The basic program format consists of a series of inquiries from 


the program concerning the values of parameters and choice of functions. 


‘The user responds to the inquiries, and, given a complete set of input 


data, the program performs the simulation and writes a report on the 
results. Tables 1 to 9 contain the results of sample sessions. The 


following discussion is organized around each of the inquiries. 
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1 Logs 

Firstly, the program edits all input information and issues 
error notifications where they are appropriate. These are all self- 
explanatory and since it is anticipated that the user will have no 
difficulty responding to them, they will not be discussed in detail. 

There is one general feature to be noted before proceeding 
with the discussion -- all inquiries which require a numerical entry, 
will require positive values. If a negative value is given (for 
example, a negative number of generations, or a negative initial 
density) the program will return to the beginning (the third inquiry). 
This permits the user to change her/his mind before entering the 
simulation. (Ordinarily, time-sharing systems provide for cancelling 
an input line before transmitting it to the system; this provision 
permits in addition, modification of previously entered material.) 
After testing for negative entries, the program also edits the input 
data for validity in,a variety of other ways. 


Program Prompts and Inquiries 


1. ***POPULATION REGULATION SIMULATOR*** 
CLEAR TABS AND SET TO 37 (RELATIVE) 
The user clears all left tabs on the terminal up to 37 spaces 
from the left margin. This is for proper alignment in the output 
table when unique control rates are entered in each generation. 


The program will inquire: 
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After each carriage return, until the user replies "YES" 


(orelY)),,. permitting, repetitive tries to set. tabs....1If they 


have not yet been set, reply "NO"(or "N''). 
SELECT A POPULATION MODEL - OS, LG, EX: 

By replying "0S" the user chooses the oscillating population 
dynamics model; "LG", the logistic, and "'EX'' the exponential. 
SELECT A CONTROL SYSTEM - IN OR SM: 

"IN" chooses an insecticidal control system; "SM", sterile 
males. 

ENTER NUMBER OF GENERATIONS: 

Reply with the number of generations for the simulation. 

This must be less than or equal to 100. A negative value returns 
to the beginning of the program, "0" is an invalid entry. 

If the insecticide control system was chosen the following prompt 
will be used: 

ENTER LD50 AND LD95 OF INSECTICIDE: 

These are to be viewed as field values; that is, the number 
of units of the material as applied in the field necessary to 
obtain 50- and 95-percent mortality, respectively. From these 
figures the program calculates the parameters of the log-dose, 
logistic-response function. 0 is an invalid entry, as is an 


LD50 greater or equal to the LD95. 
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19. 
ENTER FIXED AND VARIABLE COST OF A UNIT OF CONTROL: 

Any non-zero numbers are legal entries. Fixed cost is viewed 
as that cost of treatment which does not depend upon the dosage 
applied. Variable cost is the cost of applying a unit does 
(one sterile male in that case; one pound or ounce, or whatever, 
in the case of insecticide). The nae in the estimate of the 
user, can include any relevant values. Zero is a valid entry. 
WILL YOU APPLY UNIQUE CONTROL IN EACH GENERATION? 

Reply "YES" or "NO" depending upon whOtnet you wish to decide 
upon the control rate after seeing the projections in each 
generation ("YES"); or to specify in advance that a particular 
control rate will be used in all generations ("'NO"). 

WILL YOU USE A CONSTANT NUMBER OF STERILE MALES? 

A reply of "YES" will be followed by: 

SPECIFY NUMBER OF STERILE MALES TO BE RELEASED EACH GEN: 

Any non-negative number is valid. 

A reply of "NO" (to inquiry 9) will be followed by: 
SPECIFY PERCENT MORTALITY IN ALL GENERATIONS: 

Any non-negative number is valid. 

(A reply to "NO" to inquiry 9 will result in the program 
calculating in each generation the number of sterile males required 
to produce the percentage mortality specified in response to 
inquiry 11. A reply of "YES' will result in percentage mortality 
being calculated from the constant number of sterile males 
specified in response to inquiry 10. In each case the calculation 
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produces results which also depend upon pest density.) 
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20. 
If the sterile male control system was selected and the answer 


to inquiry 8 is "YES", the response will be: 


WILL YOU SPECIFY NUMBER OF STERILE MALES? (NOT % MORTALITY) : 


If the reply is "YES", the program will take the numerical 
input demanded during the simulation to be the number of sterile 
males released in that generation. If "NO", input will be taken 
to be percentage mortality. | 

If the insecticidal control system was selected and the same 
rate of control in each generation was specified in response to 
inquiry 8 ('"'NO''), the next inquiry will be: 

SPECIFY PERCENT MORTALITY IN ALL GENERATIONS: 

Any non-negative reply is valid. The program will calculate 
the dosage of material required to produce this level of 
mortality. Note that in the case of an insecticide, it is assumed 
that this dosage is not dependent upon density. (This is a comnon 
field assumption which has no validity. It has been shown in many 
cases that it requires far more material per unit area to obtain 
the same population control rate where density is low, than where 
it is high. However, this program does not aim at optimizing field 


application techniques). 
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16. 


17. 


zi 
SPECIFY DAMAGE FUNCTION NUMBER AND VALUE OF 100% DAMAGE: 

The damage functions are numbered as presented above in the 
texte Function numbers 1 to 4 are valid entries. The second 
entry includes all relevant costs (values) associated with 100% 
damage to the crop. 

If damage function number 1 was specified in response to 
inquiry 14: | 
ENTER PEST DENSITY AT WHICH 90% OF COMPLETE DAMAGE OCCURS: 

Any non-negative value is valid. 

If damage function number 2 was specified in response to 
inquiry 14: . 

ENTER LOG PEST DENSITY AT WHICH 90% OF COMPLETE DAMAGE OCCURS: 

Any non-negative value is valid. 

If damage function number 3 was specified in response to 
inquiry 14: 

ENTER PEST DENSITIES AT WHICH 10% AND 90% OF COMPLETE DAMAGE 
OCCURS : 

Any positive entries in which the second is greater than the 
first are valid. 

If damage function number 4 was specified in response to 


inquiry 14: 
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in the oscillating model (OS) is not normally greater than 107 


22. 
ENTER LOGS OF PEST DENSITIES AT WHICH 10% AND 90% OF COMPLETE 
DAMAGE OCCURS: : 
Any positive entries in which the second is greater than 
the first, are valid. Note that this restricts the range of 
pest density over which damage can occur to numbers greater 


than 1. 


The user will find it useful to remember that pest density 


ae 


except under unusual initial conditions and some control regimes. 
In the case of the stable models (LG and EX), the models produce 
a steady density of 100,000. The parameters of the damage 
functions can be chosen to explore a wide variety of damage 
rates, and consequently, gains and benefits from control, in 


these population density ranges. 
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21. 


ao. 
ENTER NUMBERS IN FIRST THREE GENERATIONS: 


Any positive numbers are valid. The third value is used 


_as the initial value in the stable model simulations, all three 


are used in the oscillating case. In the latter case, it will 
be useful to choose as initial values, a sequence of three from 
a simulation that has run for 20-50 years, in order to examine 
cases of regular oscillation rather tee those with a large 
transient signal effect. 

If unique control was specified in response to inquiry 8 
("YES") during execution the program will demand that the 
control percentage, or number of sterile males, be entered 
after the values of expected density and damage are written. 
The inquiry 
consists of a ''>'"' written under the column heading "CONTROL RATE". 

Following entry of the desired control rate, the carriage 
must be returned. If a negative entry is given, the program 
returns to the beginning and the current simulation is aborted. 

Following completion of the simulation summary, the user 
is asked: 

GRAPHS? 

A response of "'YES'' produces three graphs containing: 
the logarithms of uncontrolled and actual pest densities over 
time, the logarithms of expected and actual pest densities over 
time, and uncontrolled percentage damage and actual percentage 
damage over time. 


Finally, the user is asked: 
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24. 
22. FINISHED? 
A reply of nyEs" exits front the prograr, ''NO'' returns to the 
beginning for another simulation. 
Results 

The format of the results produced by the program is illustrated 
in the tables that follow. 

It is intended that questions for exploration be formulated by 
the user. Consequently, an extensive ARSE of results is not offered 
here. However, examples of the exploration of two typical questions are 
shown in Table 1 - 4, and 5 - 9. 

In the first case, the effects are shown of using sufficient 
insecticide in the oscillating population system to produce 75, 85, 
and 95 percent control in all generations (Tables 1 - 3) and the 


effects of using an undisclosed control strategy on the same system 


(Table -4). 
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Note that total short-term gains decline as increasing control 
rates are applied, whereas total long-term gains increase and then decrease, 
revealing that at some point between the 75 and 95 percent control rates 
there is an optimum level. The results (including other simulations 
on this Same system) are displayed in Fig. 9. It is clear that short- 
term gains are maximized between control rates of 50 and 85, if a control 
strategy is used which applies equal dosages of the insecticide in all 
generations. It has been shown however in Table 4 that there is some 
other strategy which produces a better result than any of those used 
above (it is left to the reader to discover this, or a better one -- 
the author would be interested in hearing about any that are better). 
It would be unusual to apply at low densities when little or no damage 
is expected. The usual procedure is to try to control under circumstances 
which produce short-term benefits -- this is left to the user to explore. 
It might be interesting to contrast the effects of controlling only at 
high density with those of applying contral at moderate, increasing 
densities and moderate, declining densities. Although the results of 
these simulations should not be taken seriously, there are some which 


‘are worth contemplating in real systems. 


Figure 9.--Short-term benefits (n) and long-term benefits (.) ina 
particular simulation (see text) as control rate in the generation is 


varied. 


Note that short- and long-term benefits are optimized at different 


control rates. 
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Control Rate 


| ieee 9,--Short-term benefits (x) and ee benefits (e) 

in a particular simulation (see text) as control rate in the 
generation is varied. | 

Note that short- and long-term benefits are Reet teat different 


control rates. 
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26. 

If the user conducts other simulations on the same system, 
allowing the value of damage, and thé pest-density damage function 
to vary, he will discover how important it is to have a detailed 
knowledge of the exact nature of the impact of the pest on the 
production of valuesin a natural resource system. The development 
of optimum forest pest management strategy requires a thorough 
integration of the results of studies in population dynamics, forest 
pest control, forest pest impact, and the dynamics of the production of 
resource values. 

Tables 5 through 8 show some interesting effects of using the 
sterile male control procedure in the exponential population system. 
It is clear that the system stabilizes at a lower pest density when 
fewer sterile males are used than are needed for extermination of the 
population and that net benefits increase remarkably with the use of 
_ large numbers of sterile males. One wonders whether there is an optimal 
trajectory towards extermination of the population using the same total 
number of sterile males. Perhaps using moderate numbers when the 
population is high and large numbers later would be better than the same 
“number in each generation, or perhaps the converse would be best. 
Table 9 shows the results of one attempt to explore this question. They 
reveal a typical result -- one can choose an optimal control procedure 
on the basis of obtaining a desired system state most rapidly, or 
on the basis of least net cost, and the two bases seldom provide the 


same result. These, however, are among topics for another discussion. 
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Copies of the program are available upon request from: 
Forest Insect and Disease Laboratory 
Northeastern Forest Experiment Station 
151 Sanford Street 


Hamden, Connecticut 06514 
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TABLE 2. 


## POPULATION REGULATIO’ SIMULATOR*#* 
CLEAR TASS AND SET TO 37 (RELATIVE) 

> 

MPD 


SELECT A POPULATION MODEL= 05,L6,EX: dos 
SOLER sAUCO) TRO lesGV ST Ee It ORS atin 
: Sg Na as (ea A ae 8 A RENEPATIONS: >29 
EheR -EDSOs a FOSS: Ore st iSEGT GUIDE speak 2 
5 ENTER FIXED AID VAR LADLE COST OF CA UNIT OF CONTROL: 91.5 
SLL YOU APPLY UNIQUE CONTROL 0H ee GEUERATION? 2 
SeeGhey PERGeriy ah WAN ye ttl AGL BRAT OMS = Digo 
SPECTFY DAMASE FUNCTION NUMBER Ae “VALUE OF 1ON™ DANAGE: >4 50 
: ENTER LOGS OF PEST DENSITIES Ar WHIGH 10% we is}toyi oll COMPLETE DANAGE OCCURS: >5 7 
ae SPECIFY NUMBERS IM FIRST THREE GENERATIONS: >500 1990 25000 : 









{ ANTHOUT CONTROL | WITH, GENT ROL: | | VALUE OF 
| | EXPECTED J CONTROL | ACTUAL | | SHORT TERH | LONG TERA fe j 
SEN DERISdTY POAMAGE | NENSITY [DAMAGE | RATE | DENSITY] DAMNAGEL DOSAGE | COST -| GAINS | SENEFAT | —-GAING 9 *SRNERO Tt lscul 
eeEKKEK ee Se ec eee Sepa ane eee eS ee ke Bee oo mnae co Cee ty ace ae 
: 1 PORE Dk: Meili O LOGS Fl eit hem gor, 0) TID? Sit My lile OD omeye car 4.88 18.37 ~49 18.37 13.49 1 
2 1N6685355. 90.5 5 7,60 55 sem G02 Smee Sore.) 101400, 10012 er 4&8 LORD 2 Toh 40.21 De aD 2 
g IZOSUGCT oe 85 09 53364. 5.8 85.00 ; 8005, Os O18) wee. 7. 4,3é 20 io 0 42.96 32.07 3 
4 GLOSTS 95 5'0:..0 2016. 0. See 8500 Baier 0.0% .777 4.83 OF Ed wo eT T 14.96 10.08 4 
cS 5 8617. al. 208, 0.0). 85. 010 as i OM een 7.7 4.86 0.01 multte: 9:1. 0.54 Sue ou: 5 
5 BSS. Dek aI 7's 9.9 85.99 Ne APL 5 Car Aze 4.38 0,01 mil ea 0.04 “4.85 6 
7 BUY, Oo1 1759, eee On 00 265, Oe OU wear Aa. 4.88 0.10 meth 5 0.01 =§5 87 7 
‘ 8 Sa 6 9 0.7 47465, 5.2 85.00 7420. 08 88 ew 77 4,88 Prorat =25 7h = az -5.00 8 
9 20d Same 210 9.5 661509. =h0. 35. 85.00 99201. De VS te i. 4,88 15% “A 10.28 on. 19 0% 3 9 
19 Saas mero on L584 Oe = 5D. Mae 855.0:0 GD 7G) Melasma 17 4.88 lise 9 14.03 34,02 205 15am ly 
11 ISHS 00. pO? 195.0 00ee 27 .Giee fion. 00 29250. Daeg ae TAT ATE 4,88 ae 221% 44,37 oo, ORF LE 
Lz LOO S12 see OD) ..0 77685, 1. One 857, 00 1168, Oe UG ieee 7 7 4,88 o.4n “448 IW 57/ 2hq Cae Le 
Ve 43742, 4.8 td 36 Oe AS aaa le] 2. OF OM ea a 7 4,88 1,92 -4.86 2.40 “2.48 13 
Yh se 0,2 159 sa « Oem oo. 00 rae OS NP stas 4,88& 0.01 -4,87 0.09 “4,80 18 
LS Sdidys 0.9 bese Daly» 85400 a0, On dle atl 4.88 0.03 “4,85 0.02 =i 870" 15 
16 LS 7 bs OZ 11555. Lok 85.00 1730. D205 Sag. 4,88 0.58 -4.50 = 0.03 “$551 _ 26 
aw B51 58. rea LLEOST 6 eke See Ga NG W1105, Od or 4,88 rel ee at 0.20 =0S 917 
1é 1933105. © Gb. 2 1329780, 956.8 85;..0:0 LO GGG: Tee E7eGiSe sit 4,88 19,54 14.66. Ley (OE LS Sdesrtie 
iN 14094620. 92.6 56574 en ee ei pag 5: (010) Bh ESh, COU AA ZATATE 4,o8 14.03 DLS 41.96 52 Vidoes 
20 GOK2759. 84.8 Bie OW «a Seem ois OD LOW7, o, G3 saa 4,86 SL =3.34 42.07 B7519 220 
De ee ee op mmbe nk ered h Aa ATERUAHE OES EADHOERSEETTATAR HD WED ANN EERADG KOREEROLERELO STA SEANAD SCRRSAA TES SRS ATS ERESAEY SETAE 
62414288. 774, 6105056.9 338. 915759.54 88.2 1535 97.65 124.7 27.06 342.9 — 245.2 


YOU PROSUCED AN AVERAGE DENSITY WHIAit WAS 1.67% CF TNAT OF TUE YNDISTURBED POPULATION, 


Ad 15.00% OF THAT OF THRE YEAR TO YEAR CONTROLLER PAPULATIO', 


YOU USED A TOTAL OF 15.53 UNITS OF NATERIAL AT A TOTAL COST OF 97.65 : 
j : THE SU NE YOUR SHORT TEAM GAINS WAS: TAN Teh TSO f\ SUeae ae i aletmay SOG 27.06 : 
THE SUI OF YOUR LONG TERM GAINS WAS: ‘342.90 FOR A NET SENEFIT OFS 245.25 
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' TABLE 3. ne 


SELECT A POPULATION MONEL= OS,LG,EX: dos 

SELECT A CONTROL eas ESO ReGh's lin 

ENTER NUMBER OF GENERATIONS: >20 

Evie DSO) ree 0 95: OF MHSECT IC ON Age 4 

ITER FIXED AND VARIABLE COST OF A UNIT OF CONTROL; >1 5 
LL YOU APPLY UNIQUE CONTROL IN EACH GENERATION? dn 


) 





9 
SPECIFY DANAGE FUNCTION NUMBER AND, VALUE OF 100% DAMAGE: >4 50 
ENTER LOGS OF PEST DENSITIES AT WHICH 10% AND 90% OF COMPLETE RAMAGE OCCURS: >5 7 
SPECIFY HOERS iN FIRST THREE GENERATIONS: >500 1000 25000 
; {| WITHOUT CONTROL || WATH CONTROL | | VALUE OF 
| | EXPECTED | CONTROL | ACTUAL | | SHORT TERM | LONG TERM 
GEN | DENSITY [OAMAGE | DENSITY [DANAGE | RATE | DENSITY] DAMAGE! DOSAGE] ~ COST | GAIN | SENEFITI] GAIN | BENEFIT [GEN 
He KH I KK KF HIKE KH HORE FHT KK TIF TE HF IRE ITE HT IT TOK ROR I IORI TOK OR RF RI HTK IK III IK RRR KKK EK HK RE KH HK ER MKRKEK KEKE RE KE 
1 LUG 6S Slo) Dl. 6 LOGS8ST.  SL6 995.00 53442. 5.16 20:00 PU 00. e227 92 abuse et 22d Lae 3d 1 
2 10664533. 90.5 L5GGo Le koe Oe to. Os GSS. = Osos 2.00 . 11.00 6.09 es SS 4G, 84 33.84 2 
3 7265467. 86.9 SLD Ls 0.6 95.00 TSIEN, 0.02 2.00% 11.00 Ong =107 81 G3.44 32.44 ES 
4 h10815. 30.0 Bie6 Oli 952100 ie 0.00 2.00 11.00 0.01 =10/99 14.98 Dea 4 
5 8317. Beek 24, lO) SHE ye 1x6, aR OBO <2 00 11.090 0,00 =~ VL. 00 0.54 -10.46 5 
55 5o5. al 108. O70 9155.00 Die 0.90 2.00 11.00 0.01 = 10% 99 0.04 -10.96 6 
if BUY, 0.1 3941. Gao, 295.00 Los 0.05 2.700 11.00 0.24 ~10.76 OlneZ =10 2918 7 
8 5169 On D7 ROU oes fe Da). 00 8570s P05 92. 00 LL. 0:0 Veod 9/6109) =0,20 =1).20 8 
9 PP SSg/ Stes PAB nl MSIF aes Se Re 59665. Go 625, 0.0 11.00 2351 I2 L292 6.98 = S02 9 
10 Sil? V2. ee sh b7 OG COS 32% 88 2195.00 23534, 25720 25.00 11.00 L502 4,02 40.67 29.67 10 
id LSOSES005: S2e2 17090. 2a Ue 95 000 855. Dea 2 92. 00 11.00 09> =J10).10i5 Leys 2) SO. e LL 
2 L8i9I6525. 65.9 331. ONG 95. 010 dis 0.90 2.00 12100 0.02 = 05918 Dlg PAS Shy 
13 43742, 4.38 2S. O20m 95.00 alr OF. 00522100 11,00 0.00 Sb St, 2.40 =. 00) eo 
14 LS a. OFZ 38. ONTO EA957.100 Die 0.00 2.00 11,00 0.00 =1 1.010 0.09 ~10.9 Leak 
15 ST 0.0 20% OF, Ve 95200 36. Ot 2.00 11.00 0.05 = 10 J9 0.02 =L0s9 8. 5 
16 L325 OC2 LOLOS. ho5 95... 010 2020., Oe caren 00 11.00 2.10 -6.90 -0.04 = Oe kG 
A 4513e. 4.9 TST Som EOL 585.00 38937. Wy. 3120 2:6 010 11.00 1987 8.87 0.51 = 10 )g 1) mene 
18 1955209. 65.2 LO 707 3 rome pe Onno O00 SoD On ts Sard 200 11.00 PAS ENS; Ldresd'S 29215 Lon Oe ke 
19 14094620. 92.6 : 87882. 8.9 95.00 4394, 0.26 25.00 11.00 4.19 =6). 31 46.01 Doe oo 
20 6042759. 84.8 LG Late Oe 2nd 5100 el. 0.01 2.00 11.00 0.10 =10.90 42.38 d2098) 20 
PHI RK EKER KARE KR ERE KEKE HRM IK KKEKKEKEMHAKEKEER erase myesastavecerarsssirestarsnretscyy $2420.55 000 ten. oi eee ane 
62414268, 774, 5045918.0 260. 252295.8 CAT Ae 40.0 220.0 125.9 -94.08 a7d01 153.1 


YOU PROMUCED AN AVERAGE DENSITY WHICH MAS 0.40% OF THAT OF THE UNDISTURSED POPULATION, 
AND 5.00% OF THAT OF THE YEAR TO YEAR CONTROLLED POPULATION, 


YOU USEN A TOTAL OF 40.00 UNITS OF MATERIAL AT A TOTAL COST OF 220,00 
THE SUN OF YOUR SHORT TERH GAINS WAS: 125.92 FOR A NET SENEFIT OF: -94.08 
THE SUI] OF YOUR LONG TERM GAINS WAS: 373.05 FOR A NET SENEFIT OF: 153.05 
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AT A POPULATION MODEL- OS,LG,EX: >os 
CT A GONTROL SYSTEM- Ut OR Siise pin 
RoNYTDER GE GENERATIONS: 220 : 
SO Lps0 S095 -0F PISECTICIDE: >. 2-2 ; 
mE EYED AMD VARIASLE COST OF iy NSE (Oe AOR WOE pe oe 

YOU APPLY UNIAVE CONTROL 0 CAG GENERATION? dy 
LEY DAMAGE FUNCTION UNSER Ald VALUE OF 100% OAHAGE: >4& 50 


oe 
ro 


LEY HUNSERS su FIRST THREE GENERATIONS: >500 1000 25000 r 


. 


62414288, 


PROPUGED AN AVERAGE DENSITY WHICH VAS 
54.58% GF THAT OF TIE YEAR TO YEAR CONTROLLED POPULATION. 

USED A TUTAL OF 1.98 UNITS OF NATERIAL AT A TOTAL COST OF 

Sui OF YOUR SHORT TERM GALS VASE 27.35 FOR A NET- BENEFIT OF: 

Su'’s OF YOUR LONG TEP GAINS WAS: 315.79 FOR A NET BENEFIT OF: 


DENSITIES AT WICH 10% AND 90% OF COMPLETE DAMAGE OCCURS: 25 7 


2S OP MAT Or ene NODLSTURSED: POPULATION, — 


12.88 
14,48 
302.92 


: - 
aus 
Ne pai aace'cseauadscensseetestbuedecntcadganedstxreesaceceueets tere tsweeeseSsee See So 0s 020 oor eee aied : 
774, 2482911.0 197. 1355157.0 142. 1.98 12.88 14.48 315.8 302.3 
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TABLE 5. 


TER HUMSER OF GENERATIONS: >10 

TER FIXED AND VARIABLE COST OF A UNIT OF CONTROL: >5 .001 

WiLL YOU APPLY UNIQUE CONTROL IN EACH GENERATION? On 

TILE YOU USE A CONSTANT MNUHSER OF STERILE TALES? dy 

SPECIFY NUNBER OF STERILE MALES TO GE RELEASED EACH GEit!: >1000 

SPECIFY DAMAGE FUNCTION NUMBER ANO VALUE QF 100% DAMAGE: >4& 500 

ENTER LOGS OF PEST DENSITIES AT WHICH 10% AND 90% OF COMPLETE DAIiAGE OCCURS: >3 5 
SPECIFY NUNSERS iN FIRST THREE GENERATIONS: >500 500 100000 





|] MATHOUT CONTROL | WITH CONTROL : | | VALUE OF 
| | EXPECTED | CONTROL | ACTUAL | ° | SHORT TERM | LONG TERM 
SEN DENSITY |OAAGE | DENSITY [DAMAGE | RATE | DENSITY! OANAGE| DOSAGE | COST) | GAIN 3] GENCEIT| GAIN 3] BENEP IT) |GEs 
ee I OR IOI III IOI II III OR IOI III IOI III IO III IIR ITOIOIOI IIIT IORI IIIT RIOT IIIT TOTTORI RI KE EERE KH 
1 100000. 90.0 1n0000. 90.0 12.96 98039. 89.83 1000. 6.00 0.86 =-5 1h 0. &6 =p 4 1 
© 2 109000. 99.0 9916. F080 1.96 97856. 89.381 1000. 6.00 0.86 ID oni 0.94 = ye Gy 2 
2 100000. 39.0 99797. 99050 196 97837. 89.81 1000. 6.00 0.86 Die at 0.95 =. U2 3 
4 100000. 930.0 Civ eag” Sil) 12.96 97855, 89.82 L000, 6.09 0.86 =D setae Oi59 Sse OS) 4 
5 109000. 90.0 Sebipkene habe 06 97834, &9,81 1000. 6.00 0.86 =5 Lh 0535 SUS 5 
6 100000. 90.90 ehelcAsseiey MAS Kui, 1.96 97854, 6981 LOGO. 6.00 1.86 =D pple 0.95 SID 6 Oe 6 
7 100000. 90.0 e997 95— a0 wo 97634. 69.81 1000. 6.00 0.86 “5/14 0.95 SAE 7 
§ 1900000. 90.0 99795. 90.0 1916 97834, $89.81 1000. 6.00 0.66 -5.14 0.95 =D «Vip 8 
9 109000. 90.0 SO 7eD 0.0 LeoG V7 Soe Os Gd LOU. 6.00 0.86 “S14 OieDi5 = 57.105 3 
10 100000. 90.90 99795. 90.0 JEG 97€34. 69.81 1000. 6.00 0.86 = Siena t 0). 915 =5.05) 190 
Fe RT RR RO IOI ORR OIRO OIRO OR ROI I ROTI TOTO RIOR TORO RIOR IR OK RIOT RIT HRT IRI IKARIA IK IRI KERR RRR HERE KEE 
1000000.0 900. 9981380.09 900. 970572.44 898, -100E 0560.00 8.597 “51.40 Ds RE, -50.62 


YOU PRODUCED AN AVERAGE DENSITY iMICH WAS 97.86% OF THAT OF THE UNDPSTURBED POPULATION, 
AO 98.04% OF THAT OF THE YEAR TO YEAR CONTROLLED POPULATION. 
YOU USED A TOTAL OF 10000.00 UNITS OF MATERIAL AT A TOTAL COST OF 60.90 


THE SUM OF YOUR SHORT TERM GAINS WAS: 8.60 FOR A HET DENERAT OF? = Ob» 4.0 


SUM OF YOUR LONG TERM GAINS MAS: 9.38 FOR A NET SENIEFIT OFF -50.62 ; 





FINISHED2>n 








oer 





ve aes 





> 
— 
| ar 

















“se 7 
: - “e ” ' a 
REM S20 - " 
| RE 30: OF KONE Foe Lem 3 2°28 FOU Y WEL CEEEYL opt -20°Cs : 2a ae 
: AOE Om OF LUnY SEL Leint yz" See Po V NEL 7EREDIL OFF 6 BEES - > ae 
eS AU DSEY ¥ LolWr or = TEED" 112 OF HVLENIYE VE W A0LVS CORL pF eo" u9 ae 
pen Wee OS" 0%s Ot SIN, OF ANE Aye 40 ABVE CURACOFTED sabat¥sion’ ) aot oA 
a Us SMORACED VA VATEVGE VENZNLA tutCN YO Ob Be? Ob Lyd ub gue NOLL ER LOMTV, 1om* 
— 7 7 4 P 7 
el - : toboode~9 200° eovTso’ oo ace* PEt oe ae tos" “T99E O2¢0700 Lars 3) -2r'eo 
f TER SORRRO enen eee set etee esnanssensenessonscevsssnweehuesotedectuns nasensennte * 
; Ts Teorag" . 20°60 $2182" 80'o tae Oas2¥" fI°ST roan" 8°90 oes = Bre 
‘ Tovpne" §30*0 3OAc2* B0'o i" te oscar" su'er yeaa" e* ou e"ee = = TP 
; 7OCdde* GO"9 easa2* d0°o fae sere’ tars? ro00* 2° on o*te Ww. 
4 245060" . ob’ o ~ Bene2” a6°O f° du di22e" go's? fano°.. e* 0 o* 97 
7 e feqapap’. 2979 92e2" B6°O IT'oe #2820* gi*ef Qee0" a) ose 
> Foote)” 86°79 aos02° se’2 f°se Sa02e° &)°8F Féve* 4 e090 o* se 
Jocoeo’  so-5 Wide" 90°06 re Oiez2" 99°47 Toue* e*o9 oO’ ee 
| 2 106898" .a9"s ation” 30'S = I°ge A427" €3°SF Taco’ . oe "se 
*$ Teoeuu" «beh gosto" 30° rae sarc’ e0°eT F600" e*00 age 
} Teegne’ jog s00000" ge°6 TaD J262u°, €2°S2 Joon’ e090 o" se _p° ee 
PE] GEABILA IDWAOE, ene Leensavetesscernenssresscntsonengnenseneoucsenseboads bisesessbonarsbenseanbenen 
; Sul = OENEILA LOWHCE) «6GENSak ovine UVSE | CENEILAl OWIVEE| voove) copL } evia | SedEElst 
: 1 FeoPtsév | euastor) youre t " | @HDML Lfen he Tow 18e 
a 1 Aen, Cow | nti Con.uar ; i ATT AE OF 


PaSCtt A AOee "a PTHZL AMUEE CEVEUVL Oe? Shug ae0 yoonDs 

#Lem TOUT Ob RETA OF IZILIED Ya MICH TOS veo 794 Gh CHMTESE Ovivee aoonys? 92 2 
PRECI PA DYUYLE FARCAION sotipey wae AVENE OF 366, BVIWeE? of 205 

GRECISA WMIZEY Ob ZAEVITS WETES 4O OF KUTEVeED BYeR esis sJou0 

“NET. AG REE V COBUAVAL HTDEM Ub Zacerire MVPeSa SA 

“FT 20M YoutA ANTONE colbisgor ti EVER eRUESVEIONL >w 

EMERY tIXED YOO AYSIWEFE COSL OF F OINEL OF COnsuur: 22 ‘Ont 

CULEG UNEDES OF CEesygtada: >fo 

SErECL ¥ COULGOF BAtLER~ fy O¢ fiw >ou 

SETECL ¥Y mOnAVi lod wowTT~ oe*re*er: pew 


TABLE 6. ; 


qa A CONTROL SYSTEIi- IN OR Silt dsm 

ENTER NUNSER UF GENERATIONS: >10 : 

ENTER FIXED ANID VARIASLE COST OF A UNIT OF CONTROL: >5 .001 

vit EL «OU pra UNIQUE CONTROL IN EACH GENERATION? On 

“VILL YOU USE A CONSTANT NUNSER OF <STERMLE HALES? dy 

SPEC LEY NUNSER UF STERILE MALES TO SE RELEASED EACH GEN: >10000 

SPECIFY OANAGE FUNCTION NUHBER AND VALUE OF 100° DAMAGE: >4 500 : 

ENTER LOGS OF PEST DENSITIES AT WHICH 10% ANC 90% OF COMPLETE DAMAGE OCCURS: >3 5 : 
SPECIFY NUMSERS (NH FIRST THREE GENERATIONS : >500 500 100000 


P HATMOUT SGONGROL 5 { WITH CONTROL | | VALUE OF P | 
| | _ EXPECTED | CONTROLI Ss sTUAL | | SHORT TERM | » LONG TERM | 
Gert DENSITY [OAMAGE] | DENSITY |OAMAGE] RATE -| DENS | DAMAGE | DOSAGE | COST | GAIN -| BENEFIT] GAtH | BENERIT 15EY 
Se ee ke og UA cee p a re PON AD EELS ED EY Do PDP EO LOPE P ERP ATE CPM rrr tere Hite Brdangeeirm cen ren 2 
1 100000. 90.0 LO000i0%" 990.107 —26;:6:7 $3550) EOS eae skOOE M05 15.00 ie 2) =i aiOiL 8.59 -6.61 1 
2 100090. 90.0 D708 50 eGo of whl nO SOU99. 18.7 a90. L005 105 15.900 §.o3 sels 274 eG Sey “4.83 2 a 
3 100099. 90.0 SG6208. nodes shi a20 AVG. gOi/ «iS LODE gOS 15.00 6.96 -6.04 10.63 Geile, bare 3 
hy 109000. 90.0 DSS90. | 8016 oe? eb 79ONSS. 87.64 E00 O05 15.00 9.00 -6.00 10.79 Cyt yant 4 
S 100000. 90.0 95908. 89.6 17.26 TI559.. 967 689 hOOR S05 15,00 9.02 = Sie 10.84 -4.1G 5 
6 100000. 90.0 SSSI, 389.6. 17.2.6 1553.6 $8 .euo) wl OOR gD 15.00 O02 = 5015 10.386 mat Ly 6 
7 100000. 90.0 SUN See AS 745 F9B25 RS Oo ok O00iE 20> 15.00 9.02 -5.98 10.86 a ee a 7 
8 100900. 90.0 Gos | od GO | wkd area ; 193522. I8/ 585 SLOOE 05 15.00 9/102 a5 1OLo 10.86 “4,14 § ; 
9 To0000'. —90'..0 Q5869.. eG lh7 nds 79322. 187 8 sOOE 205 15.00 S502 =i AG 10.86 “4,14 ) ; 
10 100009. 90.0 95369. 89.6 17.26 7Qa21. 207 a82 ehOOE S05) aaah b-20i0\ 5 9 «02 at elt 10. 86 “4.14 10 
Fe RE RIOR IO ROI III III CIOIIGIOIIOIOIOI III IOI IO IOI IOI IOI TOI ITO III IOI IIIT TTR IT IR IK RAR RIK ERK R RR KEE EERE 
1000000.0 900. 964591.44 897. 796941. 81 879. -1LOOE 06150.0 89.32 -60.68 105.1 =<hh 92 
YOU PRODUCED AN AVERAGE DENSITY WHICH WAS 79.89% OF THAT OF THE UNDISTURBED POPULATION, ° 
AID BZ Gs Or wUNl, ake THE YEAR TO YEAR CONTROLLED POPULATION. 
YOU USED A TOTAL OF ‘100000.00 UNITS OF MATERIAL AT A TOTAL COST OF 150.00 
THE SUM OF YOUR SHORT TERM GALS WAS: 39.32 FOR A MET AENRFIT OF3 60.68 ; 4 
THE SUM OF YOUR LONG TERM GAINS WAS: 105.06 FOR A HET BENEFIT OFS “44,92 
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TABLE 7. 


SELECT A POPULATION HiNDEL= 0S,LG,EX: Dex { : 
SELECT A CONTROL SYSTEIM- iN OR SM3 >sm 

ENTER NUMSER OF GENERATIONS: >10 

ENTER FIXED AMO VARIASLE COST, OF A UNIT OF CONTROL: >5 .001 

WILL YOU APPLY UNIQUE CONTROL U1) EACH] GENERATION? dn 

“ILL YOU USE A CONSTANT NUMBER OF STERILE MALES? dy 

SPECIFY NUMBER OF STERILE MALES TO 3E RELEASED EACH GEN: 425000 

SPECIFY CANAGE FUNCTION NUMBER AIO VALUE OF 100% DAMAGE: >4 500 

ENTER LOGS OF PEST DENSITIES AT WICH 10% AND 90% OF COMPLETE DANAGE OCCURS: >3 5 

SPECIFY NUMBERS iN FIRST THREE GENERATIONS: >500 500 100000 


{WITHOUT CONTROL | Witt CONTRCOL | | VALUE OF | 
| | EXPECTED | CONTROL ACTUAL | | SHORT TERM | » LONG TER | 
GENT DENSITY [DAMAGE | DENSITY |DANAGE| RATE | DENSITY| CAMAGE| DOSAGE] COST | *GAIN® | SENEFIT[ “GAIN” | DEWER«T [GEN 
TREO I ROI IO III III III IOI III COUIGIIOI ICICI OIGIGICIOIG III I IO CICIUIOIOIGICIOICIGIUITIDICI IOI ICICI ICIGIIGIIICIOICIEIEITIGI III IE IGE I I Ie 
ah 199000. 90.0 LOOOROR Ie OO C0 So.5 : CGGG75 SaaS. 250E.05 30.00 20.30 “9.70 20.30 1 al A jen 
2 100900. 99.0 SHOUNikes fale ie Siege) SEZ Lda UG eae 2 Die Oo 30.00 pee NE, = 6/05 28.42 = eno eS 2 
5 100060. 90.0 2 SSULESS So.08 50.92 Spoor Cano 0me2 J0m OS 30.00 ZG =Fs09 FF eo0 De) 3 
4 LOGOS. 9050 Sate Ci.ce mle oda Oo SPL Zee Ceo ome cO eet. 30.00 24.98 = 2 OK 37.04 7,04 4 
5 100000. 90.9 BOOCES 8729" S506" WOZG We Sx0i a2 DUicenO a 50.00 7.S\ P27) -0.78 Sao 7.05 5 
6 100090. 90.0 TOO Ste) Oo604 R7SGGe $1 26:7 4. 25018 903 30.00 Seuss 0.18 41.65 Dat lees: ) 
7 199000. 90.0 THB Wha, either Sie OY Goa Ne Cibo Gee. 7 o0eNO'S 30.00 30.94 0.94 Ue Rah! Hebrew 7 
& 109000. 90,0 TOL eee SP eb eS Vd NG OSiGre 01a take 25.0 aeO'S 30.00 SL ao Teo KG 4UG6 14.47 8 
i) 100000. 290.0 ees, EPG eh — BI) 5 HOSie C0) Nims 2 c.0emeO)o 30.90 32.04 2,04 45.48 15.48 9 
10 199000. 90.0 VEST Dee cole) Oe Oe HHE98, G0.74 .250F 05 30,00 32.44 2.44 46.31 TGeubeee LO 
Fe CE OI ROI IO EI IOI RIG ISR IOI IOISISIGOIG IIIS OGIO IGIDIGI ICICI IOI OIGIGIGIICIDIGIOIOIOIUIDIDISIIIOIIIOIIOI IID IOI OIC TOO TOT TI Ot 
1000000.0 900. 820273.87 631. SLOZ Ia b7 824, -250F 06300.0 284.9 =e OS 380.0 80.03 


YOU PRODUCED AN AVERAGE DENSITY WHICH WAS 51.02% OF THAT OF THE UNDISTURBED POPULATION, 
AN) 62.20% OF THAT OF THE YEAR TO YEAR CONTROLLED POPULATION. 





YO USED A TOTAL OF 250000.00 UNITS OF MATERIAL AT A TOTAL COST OF 300.00 
THE Suit OF YOUR SHORT TERM GAINS WAS: 284.91 FOR A NET BENEFIT OF: =115 500 
£ Sul OF YOUR LONG TERM GAINS WAS: 380.03 FOR A HET BENEFIT OF: $0.03 
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TABLE 8. eo 


ee ST A POPULATION HODEL= OS,LG,EX: >ex ; d 
ELEST A CONTROL SYSTEtI- If! OR SH: >sm ; . 

cee NUMBER OF GENERATIONS: >10 : 

ENTER FIXED AND VARIABLE COST OF A UNIT OF CONTROL: >5 .001 

WELL YOU APPLY UNIQUE CONTROL if! EACH GENERATION? >n 

“WELL YOU USE A CONSTANT MUMDER OF STERILE MALES? Dy 

SPECIFY NUNSER OF STERILE MALES .TO.SE RELEASED EACH GE: >50000 


. 


SPECIEY DANAGE FUNCTION NUMSER ANO VALUE OF 100% DANAGE: >4 500 
ENTER LOGS sé PEST DENSITIES AT WHICH 10% AND 90% OF COMPLETE DAMAGE OCCURS: >3 5 
SPECIFY NUMBERS IN FIRST THREE GENERATIONS: »500 500 100000 
ei TiOUT CONTROL | WITH CONTROL | | VALUE OF | 
| | EXPECTED |] CONTPOL | ACTUAL | | SHORT TERM | LONG TER | 
GEN | ENSITY J[OANAGE | DENSITY |OAMAGE | RATE | “SITY] DANAGE] DOSAGE] COST | GAIN | BENEFIT” GADN” | BENERIT ats 
RK EK HK KEKE KEKE KEK KEI KKH EK KEKE KREKUKERE Misarepe en reneeeres SaneinnirereaeeetrresSetrtrerrreri re retstre ttt eer trite etl c ire 
aL 100000. 909.0 foon00. 90.0 50.60 50000. 82.29 .500£°05 55.00 58.57 -16.43 3357 -10.43 2 
2 1no0nn9. 90.0 FROSicen Sot ee Oenoe : ZEOO0Se wey Oe ome OO a 55.00 55), ou OSes e 66.37 Pes Sirs 2 
3 100600. 90.0 G57 Cer eGo Oden: ANAS Wes) SOwWle Ws S50 Pip 7 20 eSe7 93.84 43.84 3 
taal 100000. 90.0 49065. 82.0 67.08 Ga SOG Vaz bee OOO 55.00 1035912 48.91 143.80 828.80 u 
5 LODOO0.~ 90.0 Se Clem Olesen ah enkah SO DG aiden OU aC 55.00 146.06 91.06 242.0% 2 15.72.04 5 
6 100000. 90.0 20705= CGu? VoZed5 Bele” ree js ions {OES 55.00 197.81 LESS Si 35 2S 9s 6 
7 100000, 90.0 8595. 46.4 92.09 680. Felis 5 0005 55.00 196.25 1G1225 bake? s a50620 7 
8 100000. 90.0 Too0% 15.5 S9seoh 28% he cirk MBO e IP 55.00 Too 20505 44YS.17 59S oLF 8 
9 100000. 90.0 70. 029 99.595 On 0.00 .500& 05 55.00 4.36 -50.64 450.00 395.00 3 
10 100000. 90.0 Oe a0 200 200 0. 0.00 .500E 05 55.00 0.01 -54.99 450.00 3957010 5 #10 
HK KE KH KE KKK KKK RR KK KH KEKE KKK REE HERI EEK IKKE KK EHR EKER KR EK KARE KE OI OI OR ROR TOR OR IO IORI ROI OR IRR ROTOR TI RIO OR IORI RIOR OR eR 
1000090.0 900. 357379.94 S523 138893.12 Biel ~500E 06550.0 893.5 Sipois 2637. 2087. 


YOU PRODUCED AH AVERAGE DENSITY WHICH MAS 13.89% OF THAT OF THE UNDISTURBED POPULATION, : 
A‘lD 38.86% OF THAT OF THE YEAR TO YEAR CONTROLLED POPULATION. 

YOU USED A TOTAL OF 500000.00 UNITS OF NATERIAL AT A TOTAL COST OF 8 550.00 

THE SU! OF YOUR SHORT TERN GAINS WAS3 893.52 FOR A NET BENEFIT OF: Soe 

TUE SUN. OF YOUR LONG TERM GAINS MAS: 2657 .0de COR AM hed OSENERIT O's 2087.01 
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: TABLE 9. 






3 : 
SELECT A POPULATION MOMEL= 0S,LG,EX: ~ 
ELECT A CONTROL svedene TORE Shs e2s 
ENTER NUMNSER OF GENERATIONS: >10 
ENTER FIXES AND VARIABLE COST OF A UNIT OF CONTROL: >5 .001 
Wilk YOU APPLY UNIQUE CONTROL IM EACH GENERATION? >y 
WiLL YOU SPECIFY NUNSERS OF STERILE HWALES? (OT % MORTALITY): Oy 
‘ SPECIFY DANAGE FU: eTION HUMBER’ ANO VALUE OF 100% DAIIAGE: >4 500 
ENTER LOGS OF PEST DENSITIES AT WHICH 10% AND 90% CF COMPLETE DAMAGE OCCURS: 93 5 
SPECIFY NUMSERS IN FIRST THREE GENERATIONS: >500 500 100000 8% 
| WITHOUT CONTROL | WITH CONTROL | | VALUE OF | 
| | EXPECTED | CONTROL} ACTUAL | [o> SHORT TE Re | LONG TEN 
GEN | DENSITY | DAMAGE | DENSITY |OANAGE RATE | DENSITY] GANAGE| DOSAGE | COST AEGAN Le GENE RIN (mG AIH tee | mecieiy Eesti [pee 
AE Pdi ee OA ee asda ie ages Ant ne IA OO AE Ee TON PREAH ee POO REE EET OP Yr reer Set TE y oee CL 
1 100000. 90.0 100000. 90.0 >100000 ‘ s 
66.67 BoE, Tein ley oats 195.00 70.35 SMS LE) USES -34.65 1 
2 100000. 90.0 61401. 85.0 >200000 
66.69 Slice US. 205. 200k 0G 205.00 198.84 -6.16 224.02 197.02 2 
3 1090000. 510.10 - 18954. 64.8 >1509000 ° 
94.06 PEZGse ke deer DOE NOG 155.060 268.64 113.64 394.66 233.66 5 
4 100000. 90.0 PE PAPA PEG NUTIG 
OS a0 124. T3o9, S00E 05 poo. 00 LOS) 71.59 442.55 407.55 4 
; 5 1000090. 90.0 309. 3.5 >10000 
98.48— De 0.07 .100E 05 15.00 17.14 2.14 449.67 434.07 5 
is) 100000. 90.0 12s 0.2 5000 ; 
99.88 0. 0.00 .500£ O04 10.00 0.80 =O) 20 450.00 440.00 6 
7 100000. 90.0 0. 0.0 >3090 
100,00 0. 0.00 .300F O04 8.00 0.00 -8.00 450.00 442,90 7 
$ 100000. 90.0 De O07 0 
00 : 0. 0.00 .0 0.0 0.0 0.0 450.00 450.90 8 
3 100000. 90.0 0. 0.9 >2000 


POPULATION EXTERMIMATED 


Fe ccc gud ahiGdenderWGEUEndPEstbTRANESRERADECETALAGHERTRAREOREDERS EROS NESTERRS GEER RAS ORE RENE CREDSET EEO A SENS AASR CSS EAS 


800000.00 720. 183461.44 266. 42756.875 13h. -498E 06533.0 662.4 , 129.4 295 da 2398. 





YOU PPONDUCED A! AVERAGE DENSITY WHICH WAS 5.34% OF THAT OF THE UNDISTURSED POPULATION, 
AD 23.31% OF THAT OF THE YEAR TO YEAR CONTROLLED POPULATION, 
YOU USED A TOTAL OCF 498000.00 UNITS OF MATERIAL AT A TOTAL COST OF 533.00 
TYE SUI OF YOUR SHORT TERM GAINS WAS: H62,090 .-0R AT HET SEVERIT OFS 129.36 
= Tii— SUI OF YOUR LONG TERM GAINS. WAS: 2931.24 FOR A SIET GENEFIT OF: 2398.24 : 
FINISHED?>n : 
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